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General  Introduction 

In  the  period  since  submission  of  the  last  report  the  grant 
sponsor  has  not  required  any  analyses  at  Salford  of  materials 
developed  at  Wright  Patterson  AFB.  Consequently,  the  period 
has  been  employed  to  conduct  further  back-up  research  relevant 
to  the  needs  of  the  sponsors  in-house  programme.  This  study 
divides  into  four  areas  as  follows: 

1.  Analytical  analysis  of  the  data  derived  from  the  use 

of  Rutherford  backscatteri no  -  channelling  measurements 
of  disorder  production  in  ion  implanted  semiconductors 
for  the  determination  of  damage  production  mechanisms. 

2.  The  redi stri buti on  and  lattice  i ncorporati on  of 
implanted  impurities  in  Si  during  controlled  furnace 
annealing  to  reorder  the  implantation  damaged  Si. 

3.  The  production  of  disorder  in  InP  by  light  and  heavy 
ion  implantation  as  a  function  of  implant  flux  and 
fluence  and  implant  temperature. 

4.  Investigation,  in  parallel  with  3),  of  the  damage 
creation  and  annealing  processes  associated  with  the 
interaction  of  the  RBS/channel 1 i ng  analysis  ion  probe 
with  heavy  ion  implantation  damage  in  InP. 


Substantial  progress  has  been  made  in  all  four  areas  of  study. 


which  are  discussed  separately  and  sequentially  in  the  following 
report,  and  investigations  in  areas  1)  and  2)  have  reached  a 
satisfactory  conclusion.  Work  in  areas  3)  and  4)  is  continuing 
at  Salford,  under  alternative  funding  arrangements  since  USAF 
support  is  now  terminated,  because  of  the  importance  of  under¬ 
standing  ion  implantation  processes  in  InP  and  related  materials 
which  have  great  potential  application  in  opto-el ectroni c  devices 


THE  ANALYSIS  OF  RUTHERFORD  SCATTERING-CHANNELLING  MEASUREMENTS  OF  DISORDER 
PRODUCTION  AND  ANNEALING  IN  ION  IRRADIATED  SEMICONDUCTORS 


Rutherford  scattering  and  channelling  of  light  probe  ions  (e.g.  He+)  has  been 
extensively  used  for  studies  of  disorder  production  in  ion  implanted  semiconductors 
The  technique  is  employed  to  determine  the  number  of  lattice  atoms  displaced 
more  than  ^0.2  8  from  equilibrium  positions  but  unlike  microscopic  techniques  does 
not  give  direct  information  about  the  nature  of  defects  (although  inferences  can  be 
made  from  dechannelling  studies^),  such  as  if  the  defects  are  isolated  or  in  agg¬ 
lomeration  as  in  amorphous  zones.  Study  of  the  observed  backscattering  signal, 
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either  from  fixed  depth  or  depth  integrated'  ’  as  a  function  of  implant  ion 

fluence,  has  however  been  used  to  infer,  qualitatively,  the  nature  of  the  disorder 

(2) 

production  process  for  a  variety  of  energy  deposition  density  conditions  in  Si'  ' 

and  GaAs^  implanted  at  ^40  K  to  minimise  annealing  processes.  These  processes 

range  from  simple  defect  accumulation  until  a  sufficiently  high  local  defect  density 

is  reached  with  increasing  fluence  that  relaxation  to  amorphousness  occurs^’^  to 

direct  impact  amorphisation  in  individual  cascades  where  defect  densities  are 

(81 

sufficiently  large  to  cause  instant  collapse  to  amorphousness'  ' .  Various  authors 
have  analysed  these^  and  composite  extended  models^’^’^  of  amorphousness 
accumulation  and  Carter  and  Webb^0’^  have  indicated  the  general  difficulties  in 
assessing  disorder  production  models  from  RBS/channell ing  studies  if  the  production 
modes  are  complex  and  the  manner  in  which  the  technique  responds  to  different  defect 
structures  is  unspecified.  For  less  complex  disorder  production  modes  and  by  making 
reasonable  assumptions  about  the  technique  response  however,  some  insight  into  the 
form  of  backscattering  yield  -  ion  implant  fluence  functions  can  be  obtained  as  is 
discussed  in  the  present  communication.  It  thus  becomes  possible  to  infer  the 
importance  of  different  disorder  generation  processes  from  RBS/channelling  -  ion 


fluence  studies.  It  will  also  be  shown  how  simple  annealing  processes  modify  dis¬ 
order  accumulation  and  thus  again  how  the  operation  of  such  processes  may  be  inferred 
from  RBS/channel 1 ing  -  ion  fluence  measurements. 


For  simplification  and  definiteness  we  assume  that  an  implant  ion  may  produce  a 
completely  amorphised  area  of  solid  (at  depth  z),  and  an  area  ad  which  is  dis¬ 
ordered  but  not  amorphised.  If  in  further  ion  impact  parts  of  the  areas  afl  or  ad 
produce  further  disorder  in  an  already  disordered  but  not  amorphised  area  at  depth  z, 
the  area  of  overlap  is  transformed  to  amorphousness.  This  is  the  single  overlap 
model  of  earlier  authors^*^’^*^  and  may  readily  be  adapted  to  multiple  overlap 
situations  and  v/here  ad  is  non  discrete  but  composed  of  areas  of  variable  disorder 
(10’1D.  A  set  of  coupled  differential  equations  describing  the  fractions  of  area  at 
depth  z,  amorphised  Aa>  disordered  Ad  and  ordered  Au  is  readily  developed  and  solved 
(7,9,10,11)  as  a  func-t i on  of  ion  fluence 


Thus, 


=  1  -  (1  + 

ad$)  exp  -  (aa  +  ad)$ 

(la) 

Ad 

=  ad$  exp  - 

(aa  +  ad)$ 

(lb) 

and 

\ 

=  exp  -  (aa 

+  ad>* 

(lc) 

where 

Aa 

+  Ad  +  Au  ' 

1 

(Id) 

In  RBS/channel 1 ing  studies,  the  backscattering  yield  from  depth  z  is  composed 
of  projectiles  passing  through  the  solid  in  a  random  direction  scattering  from  all 
atoms  and  projectiles  executing  an  aligned  trajectory  scattering  from  atoms  displaced 
by  >0.28  Thus  the  yield  Yq  from  the  partially  amorphous,  partially  disordered 

and  partially  ordered  crystal  at  depth  z,  compared  to  the  yield  from  a  totally 
amorphised  solid  is  given  by 

VVR  ■  Aa  +  Ad  0  -*>>r+  V  V  <2> 

S 

where  x  is  the  random  fraction  of  the  analysis  beam  at  depth  z,  Nd  is  the  disorder 
density  in  disordered  crystal  and  N  the  target  density. 


In  equation  (2)  the  atomic  density  of  amorphised  areas  is  assumed  equal  to  that 


-  - 

•  *  ■  * 
k*».' 


Bit 


il 


of  the  undamaged  substrate.  This  is  of  course  slightly  incorrect  and  will  lead  to 

multiplication  of  the  first  term  on  the  right  hand  side  of  equation  (2)  by  a  constant 

slightly  less  than  unity.  If  N.  is  the  critical  local  defect  density  to  result  in 

Ndc 

collapse  to  amorphousn’ess  then  we  may  write  ^ —  =  f  as  the  critical  fractional  defect 

density  to  lead  to  collapse.  Further  if  single  overlap  conditions  are  assumed  then 
Nd 

Nd  =  ar,d  ec!uat‘’on  (2)  may  be  rearranged  using  equation  (Id),  to  give, 

Yn 

y  =  (1  -  x)Aa  +  (1  -  x)f/2  Ad  +  x  (3) 

R 

Substitution  of  equations  (la  -  lc)  into  equation  (3)  leads  to  the  result: 


Mi 


K 

-  A  -■ 


v- 
•*>. 
<\v.  . . 
v  •  ■ 


>3 


Yn  f 

y-  *  1  -  (l-x)exp  -  (a.+ad)*  +  <(l-x)(i-l)  ad«  exp  -  (aa+ad)») 

R  (4) 

Since  the  yield/fluence  function  is  usually  determined  over  several  orders  of  magni¬ 
tude  in  fluence,  $,  it  is  usual  to  plot^’2^  the  log  YD/YR/log$  function  and  observe 

(2  9  10) 

any  features  of  the  slope  of  this  function v  *  *  ' .  It  is  readily  deduced  from 

d(log  Yd/Yr) 

equation  (4)  that  this  slope  — d~  $ - =  m. 


{(l-x)(aa+ad)exp-(aa+ad)<M(l-x)(f/2-l)Had-ad(aa+ad)4>}exp-(aa+ad)$} 


tl-(l-x)exf 


For  small  fluence,  $,  this  equation  may  be  expanded  to 


w} 


m  =  <J> 


{(1-x)(aa+ad)+kad‘kad(aa+ad^} 


-x  +(Vad)9+kad^  +FIaa+ad, 


•  M 


where  k  =  (1-x) (f/2-1 ). 

For  small  fluence,  and  consequently  small  disorder  levels,  X<<1  and  equation  (6) 
becomes 

(a  +a.)  +  2kaH  1 

m  *  1  -  *  <7) 

As  fluence,  $,  tends  to  zero,  the  slope  m  tends  to  unity,  typical  of  all  disorder 
models  and  employing  RBS/channel  1  ing  analysis^*^ . 


We  note  some  special  cases  of  equation  (7).  First  if  direct  impact  amorphisa- 
tion  is  the  dominant  disordering  process  than  aa>>acj»  and  as  fluence  increases,  the 
slope  always  decreases  from  its  initial  value  of  unity.  Such  behaviour  is  observed 


with  high  energy  density  cascade  processes 


(2,5) 


Secondly,  if  disorder  overlap 


processes  are  dominant  then  a(j>> a&  and  since  k  =  f/2  -  1,  then  if  f  is  a  small 

fraction,  k  may  be  negative  and  the  slope  increases  with  fluence,  $.  This  behaviour 

(2  51 

is  observed  with  low  energy  density  cascade  processesv  ’  At  higher  fluences  the 
approximation  of  equation  (7)  becomes  invalid  and  in  all  cases  m  finally  tends  to 
zero  as  *-*■*.  Thus  low  energy  density  cascades  can  lead  to  sigmoidal  behaviour  of 
the  log  Yp/Yp/log'J’  function. 

It  is  also  interesting  to  note  that  departures  from  slope  linearity  may  be 
expected  to  occur,  from  equation  (7),  when  . 


(a  +a.)  +  2ka. 

d  .  L.-  <>.«„)  *  2 


ivv +  kad  3  0 

For  direct  impact  amorphisation  dominance,  a  >>a  .,  this  indicates  4>aa^2.  If  disorder 

accumulated  linearly,  without  overlap  of  already  amorphised  regions,  then  #aa^l 

would  be  the  saturation  amorphous  level,  so  that  one  concludes  that  linear  slope  is 

preserved  up  to  a  substantial  fraction  of  complete  amorphisation,  as  noted  for  high 

energy  density  cascades ^ .  If  disorder  overlap  dominates,  acj>>aa»  then  ^ 

and  again  assuming  linear  accumulation  as  a  guide,  this  indicates  that  departure  from 

YD 

linearity  occurs  when  the  measured  disorder  level  y—  is  of  order  f/2.  Thus  if 

R  (31 

collapse  to  amorphousness  occurs  when  the  local  defect  density  is  of  order  5-10%  v  ', 

the  backscattering  yield/f luence  function  departs,  on  double  logarithmic  plot,  from 

linearity  when  the  measured  disorder  is  at  about  this  level.  Thus  the  inference  of 

Ref.  2  that  superlinear  behaviour  of  the  log  Yp/Y^/log'f  fluence  function  at  about  the 

10%  disorder  level  indicates  a  dominance  of  simpler  defect  production  with  subsequent 

amorphous  collapse  is  shown  to  be  valid. 


Where  the  relative  values  of  afi  and  a^  are  less  well  specified  or  more  complex 


models  of  disordering  are  assumed,  the  resulting  predictions  of  the  behaviour  of  the 
slope  m,  although  analytically  tractable  lead  to  cumbersome  results  which  are  difficult 
to  compare  with  experiment^0’^ ^ .  The  present  treatment  for  more  clear  cut  cases 
does  suggest,  however,  that  useful  indications  of  dominant  disordering  modes  may  be 
inferred  from  slope  changes,  and  the  disorder  level  at  which  they  occur. 

Finally,  we  note  the  complications  which  will  arise  from  concurrent  annealing 
processes,  such  as  will  occur  during  most  non  low  temperature  implantation  of  semi¬ 
conductors^.  For  brevity  we  will  specify  the  form  of  disorder  production  as  direct 
impact  amorphisation  alone  and  consider  a  first  order  annealing  process  described  by 
a  single  time  constant  t  (which  will  be  temperature  dependent).  The  rate  equation  for 
amorphousness  production  is  then 

dAa  Aa 

_j|  =  Jaa  <1-Aa)  -  JL  (9a) 

where  J  is  the  ion  flux  density. 

This  equation  may  be  rewritten  in  terms  of  ion  fluence  $  as 


=  aa 


11  -  Aa  O  +  iT7>' 


which  solves  to 

a,JT 

A.  = 


1  +  a  Jt 

J  {i  .  exp  -  (-  -)aa*3 


The  RES/channell ing  yield  will  be  approximately  proportional  to  this  area,  and  so  it 
is  readily  deduced  that  the  low  fluence  slope,  m,  of  the  log  Yp/Yp/log^  function  is 

%  1  -  W1  <’°) 

Just  as  for  the  previous  cases  considered,  the  very  low  fluence  slope  is  unity 
and,  as  for  the  non-annealing,  direct  impact  amorphisation  case,  the  slope  m  decreases 
with  increasing  fluence  The  faster  the  annealing  process  (smaller  t)  the  lower  the 
fluence  (and  associated  disorder  level)  at  which  reduction  below  unity  slope  occurs. 

The  larger  the  ion  flux  density,  the  larger  the  disorder  level  at  which  reduction 


EyyniTVUm  *  FfTWfTW  U.'«  H ■  vr  VV  %  -  rv  ,r.  -  -■  w  -  -  ■ 


below  unity  slope  occurs.  It  is  thus  partly  clarified  why  comparison^ "of  Ar  ion 
implantation  of  Si  at  40  K  and  at  300  K  lead  to  a  reduced  log  disorder  -  log  fluence 
slope  at  the  higher  temperature  at  all  fluences.  In  this  case  however  some  disorder 
accumulation  process  probably  occurs  also,  so  that  the  full  analysis  is  rather  more 
complex  than  the  above.-  In  such  cases  the  details  of  the  disorder-fluence  behaviour 

( 13) 

depend  sensitively  on  the  model  assumptions v  '  and  if  both  disordering  and  annealing 
processes  are  complex,  so  also  is  the  behaviour  of  the  slope  m.  If,  however,  low 
temperature  studies  suggest  dominance  of  direct  impact  amorphisation,  then  higher 
temperature  studies  and  variable  flux  density  studies  will  allow  estimation  of  the 
annealing  time  constant,  t}  if  single  valued,  from  observation  of  the  stages  at  which 
m  departs  from  unity  for  different  flux  densities. 

It  should  also  be  noted  that,  as  a  further  indication  of  annealing  processes,  the 
saturation  disorder  level,  as  predicted  from  equation  (9c)  is  an  increasing  function 
of  flux  density  and  a  decreasing  function  of  temperature  (decreasing  t). 

It  is  thus  clear  that  study  of  RBS/channell ing  results  of  log  (yield)  as  a 
function  of  log  (implant  fluence)  can  be  used  to  distinguish  probable  disordering 
mechanisms,  to  estimate  disorder  densities  required  for  amorphousness  collapse  and  to 
investigate  thermal  annealing  processes.  As  noted  earlier  more  complex  model  assumptions^ 
(10,11,13)  can  ^  treated  similarly  but  with  less  clarity  of  prediction. 
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EPITAXIAL  RECRYSTALLISATION  OF  GALLIUM  IMPLANTED  (100)  SILICON 


ABSTRACT 

High  depth  resolution  Rutherford  backscattering/ channelling  and  low  angle 
X-ray  texture  camera  analysis  have  been  employed  to  study  the  recrystall¬ 
isation  behaviour  of  gallium  implanted  (100)  silicon  during  low  temperatur 
furnace  annealing  at  465°C. 

The  recrystallisation  behaviour  was  found  to  depend  strongly  on  the 
gallium  concentration,  exhibiting  three  distinct  recrystallisation  regimes 

For  low  fluence  gallium  implants  peak  concentration  ''’0.3  atomic  percent, 
epitaxial  recrystallisation  was  observed  to  proceed  to  completion  with  a 
substantial  fraction  of  the  implanted  gallium  being  incorporated  on  to 
substitutional  lattice  sites.  The  epitaxial  recrystallisation  rate  was 
enhanced  by  gallium  concentrations  ^0.2  atomic  percent  and  substitutional 
concentrations  in  excess  of  the  maximum  equilibrium  value  were  observed. 

In  addition.  X-ray  analysis  implies  the  existence  of  a  thin,  <3  nm, 
preferentially  oriented  polycrystalline  surface  layer  after  annealing. 

For  higher  gallium  fluences,  peak  concentration  '''1.8  atomic  percent 
epitaxial  recrystallisation  no  longer  proceeded  to  completion  but  ceased 
some  20  nm  from  the  silicon  surface.  Considerable  gallium  was  observed 
to  be  redistributed  towards  the  surface  by  the  advancing  crystal-amorphous 
interface  and  the  epitaxial  recrystallisation  rate  was  again  observed  to 
be  enhanced  by  gallium  concentrations  <0.2  atomic  percent,  reaching  a 
maximum  value  of  7.5  times  that  of  undoped  amorphous  silicon  layers. 
Substitution  gallium  concentrations  v3.o  x  10^°  Ga/cm^  were  observed,  a 


value  comparable  to  that  measured  following  pulsed  laser  annealing 


Finally,  X-ray  analysis  confirms  that  the  20  nm  surface  region  contains 
preferentially  aligned  polycrystallites. 

For  the  highest  gallium  fluences  investigated,  peak  concentration  M.O 
atomic  percent  epitaxial  recovery  was  initially  slow  and  extensive  gallium 
redistribution  was  observed  after  ''-IS  minutes  annealing  X-ray  analysis 
again  confirmed  the  presence  of  preferentially  oriented  polycrystallites. 


INTRODUCTION 


A  chin  amorphous  silicon  layer  can  be  produced  on  a  single  crystal 

silicon  substrate  by  high  fluence  ion  irradiation.  Such  layers  are 

mecastable  and  recrystallise  epitaxially  from  the  underlying  substrate 

when  heated  to  temperatures  around  450°C*.  The  epitaxial  regrowth 

kinetics  of  uncontaminated  layers  have  been  investigated  in  detail  and 

the  recrystallisation  process  is  observed  co  be  thermally  activated  with 

12  3. 

a  well  defined  activation  energy  ’  ’ 


The  activation  energy  is  observed  to  be  a  constant  for  all  crystal 
orientations  but  the  regrowth  velocity  is  orientation  dependent,  being 

twenty-five  times  faster  in  the  <100>  direction  than  in  -he  <lil>  airection 

o  4  .  , 

at  530  C  '  Thermodynamic  arguments  and  modelling  of  tne  amorphous-crystal 

interface  have  led  to  a  reasonable  understanding  of  the  recrystallisation 

process  and  its  orientation  dependence^'  The  recrystallisation  interface 

is  believed  to  break  up  into  terraced  (111)  planes,  the  number  of  terraces 

being  determined  by  the  angle  the  regrowth  plane  subtends  to  the  (111)  surface 

Epitaxial  recrystallisation  is  then  thought  to  proceed  via  the  growth  of 

<110>  ledges  on  this  terraced  interface.  The  orientation  dependent  regrowth 

velocity  then  simply  reflects  the  ledge  concentration  (or  more  particularly 

the  kink  concentration) . 


In  addition  to  orientation  dependence  the  epitaxial  recrystallisation  velocity 
of  silicon  displays  a  strong  dependence  on  the  type  and  concentration  of 
impurities  present  during  regrowth.  Low  concentrations  (y  0.5  atomic  perceiT 
of  certain  impurities  are  observed  to  enhance  or  retard  the  epitaxial  regTov-A 
velocity  by  up  to  an  order  of  magnitude.  Initial  investigations  by  Csepregi 
et  al"’''  suggested  that  electrically  active  (B,P,As)  impurities  were 
responsible  for  eiihancing  epitaxy  while  electrically  inactive  species  (N,0,C) 


retarded  it. 


More  recent  data  by  Suni  et  al  and  Lietoila  et  al  supports  the  view  that 

it  is  the  electronic  properties  of  the  impurity  which  determines  the 

rate  of  epitaxy.  A  model  for  this  behaviour  has  recently  been  proposed 

9 .  - 

by  Williams  and  Elliman 

For  high  impurity  concentrations  C>  solubility  limit)  epitaxial  regrowth 
is  considerably  more  complex  and  both  electrically  active  and  inactive 
species  are  observed  to  retard  epitaxy*® ’ ** ’  In  cases  where  epitaxy  is  able 
to  proceed  impurity  redistribution  often  results  as  a  consequence  of 
segregation  at  the  advancing  amorphous-crystal  interface,  and  in  extreme 
cases  epitaxy  completely  ceases  and  polycrystalite  nucleation  may  ensue**’ 

To  understand  the  role  of  impurities  in  the  recrystallisation  kinetics  of 
silicon  it  is  essential  to  establish  a  data  base  of  information  from  which 
comparisons  and  patterns  can  be  explored.  At  present  detailed  re- 
crystallisation  kinetics  exists  only  for  a  few  species  and  consequently  it 
is  difficult  to  correlate  observations  and  theoretical  predictions.  This 
paper  presents  the  results  of  a  detailed  investigation  into  the  recrystall¬ 
isation  behaviour  of  gallium  implanted  (100)  silicon.  Regrowth  kinetics 
and  impurity  redistribution  are  examined  in  decail  as  a  function  of  gallium 
concentration  in  the  range  up  to  4.0  atomic  percent, 

EXPERIMENTAL 

High  depth  resolution  Rutherford  backscattering  (RBS)  and  channelling  of 
2  MeV  He+  was  used  in  conjunction  with  X-ray  texture  camera  analysis  to 


characterise  the  recrystallisation  behaviour  of  gallium  implanted  (100) 
silicon  during  annealing  at  465°C 


Commercially  prepared,  float  zone  (10O)  silicon  substrates  with  low  bulk 
impurity  concentrations  (>100  k ft  cm)  were  employed  throughout  this  investi¬ 
gation.  All  samples  were  chemically  cleaned  and  the  native  oxide  removed 
just  prior  to  implantation.  During  implantation  the  target  chamber 
pressure  was  maintained  at  ''-10  ^  torr  and  samples  were  cooled  to  <150  K  by 
a  liquid  nitrogen  cold  finger.  All  samples  were  tilted  7?  from  the  incident 
beam  direction  to  minimize  ion  channelling  effects  and  all  implants  were 
performed  with  40  kV  Ga+  ions. 


Post  implant  annealing  was  carried  out  in  a  quartz  tube  furnace  with  a 
flowing  atmosphere  of  85%  and  15%  H^.  The  furnace  temperature  was 
profiled  with  a  Ni/Cr  :  Ni/Al  thermocouple  and  samples  were  always  positioned 
within  the  accurately  determined  constant  temperature  zone.  Minimum  anneal 
times  of  lO  minutes  were  employed  to  reduce  the  anneal  time  uncertainty  due 
to  warm  up  and  cool  down  periods. 


RBS  and  channelling  measurements  were  performed  with  2  MeV  He+  and  a  solid 

state  detector  (15  keV  fwhm)  was  employed  in  a  low  angle  exit  geometry  to 

provide  an  absolute  depth  resolution  of  ^80%  ^ .  The  actual  scattering 

geometry  employed  is  depicted  inset  in  subsequent  figures.  RBS  depth 

scales  were  determined  from  the  semi-empirical  stopping  powers  of  Ziegler 
18 

and  Chu  and  gallium  concentrations  were  determined  by  assuming  a  constant 

22  -3 

atomic  density  of  5  x  10  cm  for  silicon. 


X-ray  analysis  was  performed  with  a  low  angle  cylindrical  texture  camera  ' 
employing  Cu-KO  radiation.  The  collimated  X-ray  beam  was  incident  at  12° 
to  the  sample  surface  to  sample  shallow  surface  layers.  The  geometry  of  the 
camera  was  such  that  the  X-ray  beam  was  coaxial  with  the  cylinder  axis  and 


the  sample  was  positioned  at  the  centre  of  the  cylinder  on  a  rotating  shaft. 


RESULTS  AND  DISCUSSION 


As  outlined  above  the  epitaxial  recrystallisation  process  is  well 
character isted  for  uncontaminated  amorphous  silicon  layers  and  a 
reasonable  understanding  of  the  recrystallisation  mechanism  is  emerging. 

For  contaminated  layers  however,  the  available  experimental  data  is  limited 
and  consequently  the  physical  understanding  of  impurity  effects  is  less 
complete. 

The  present  discussion  examines  the  recrystallisation  kinetics  of  gallium 

•  O 

implanted  (100)  silicon  during  thermal  annealing  at  465  C.  The  results 
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for  low  (6  x  10  Ga/cm  ),  medium  (3.2.  x  10  Ga/cm  )  and  high 
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(7.5  x  10  Ga/cm  )  fluence  gallium  implants  are  summarised  in  figures 
1,2  and  3  respectively.  The  top  half  of  each  figure  depicts  the  as 
implanted  (full  curve)  and  post  anneal  (broken  curve)  gallium  distributions 
with  the  final  substitutional"  gallium  distribution  indicated  by  the  shaded 
area  of  the  curve.  The  lower  portion  of  each  figure  depicts  the  regrowth 
process  of  the  amorphous-crystal  interface  for  isothermal  annealing  at  465°C. 
It  should  be  noted  that  the  terms  low,  medium  and  high  fluence  referred  to 
above  are  arbitrary.  They  have  only  been  employed  to  simplify  the  following 
discussion  and  refer  to  the  implants  reported  in  figures  1-3. 

For  the  low  dose  implant,  figure  1,  epitaxial  regrowth  is  observed  to  re¬ 
distribute  a  small  fraction  (y  13Z)  of  the  implanted  gallium  towards  the 
silicon  surface.  The  peak  gallium  concentration  is  reduced  from  'v  0.3  to 
0.25  atomic  percent  and  the  final  distribution  is  observed  to  have  a  narrow 
gallium  peak  at  or  near  the  silicon  surface.  The  width  of  this  peak  is 
within  the  RES  system  resolution  and  its  exact  height  and  depth  distribution 
cannot  be  extracted  from  the  data.  Observation  of  the  gallium  distribution 
.as  a  function  of  anneal  time  shows  that  gallium  is  redistributed  by  the 
advancing  recrystallisation  front,  as  previously  reported  for  indium  implanted 


silicon 


This  form  of  redistribution  is  commonly  observed  for  impurity 


11. 


concentrations  above  the  equilibrium  solubility  limit  and  it  has  been 
speculated  that  interfacial  stress  is  the  driving  force  for  the  re¬ 
distribution^  '  During  regrovth  a  high  fraction  of  the  implanted  gallium 

is  incorporated  on  to  substitutional  lattice  sites  and  a  peak  substitutional 
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gallium  concentration  of  'v  1  x  10  Ga/cm  is  measured  after  regrowth. 

This  value  is  v  2.5  times  greater  than  the  maximum  equilibrium  substitutional 
solubility  of  gallium  in  silicon  (4  x  10^  Ga/cm^)^" 

From  the  lower  half  of  figure  1  the  average  epitaxial  regrowth  rate  of 
silicon  is  observed  to  be  enhanced  by  the  presence  of  low  gallium  con¬ 
centrations.  The  average  regrowth  rate  for  the  case  shown  in  figure  1  is 
0.28  nm/min,  compared  to  v  0.20  nm/min  for  undoped  amorphous  silicon. 

The  regrowth  rate  is  clearly  concentration  dependent,  reaching  its  peak  value 
of  0.55  nm/min  at  the  peak  of  the  gallium  distribution.  For  gallium 
concentrations  below  -v  0.2  atomic  percent  the  regrowth  rate  remains  essentially 
at  the  undoped  level  of  'v  0.2  nm/min  while  for  concentrations  above  this 
value  the  rate  increases  significantly,  remaining  approximately  constant 
throughout  the  broad  peak  in  the  gallium  distribution.  The  regrowth  rate 
finally  slows  again  as  the  recrystallisation  front  passes  through  the  gallium 
peak  and  approaches  the  silicon  surface.  From  figure  1  it  is  clear  that 
epitaxial  recrystallisation  approaches  completion  after  annealing  at 

O 

465  C  for  180  minutes.  Final  Rutherford  backscattering  analysis  suggests  that 
the  recrystallised  layer  is  near  perfect  crystal  as  shown  in  figure  4. 

This  matter  will  be  discussed  later  with  regard  to  X-ray  texture  camera 
analysis.  It  should  be  noted  that  the  maximum  regrowth  rate  of  0.55  nm/min 
observed  above  is  approximately  three  times  that  observed  for  undoped  amorphous 
layers,  an  enhancement  which  is  similar  to  that  observed  by  Csepregi  et  al^ 
for  0.22  atomic  percent  phosphorus  impurity  during  annealing  at  475°C. 


The  regrowth  behaviour  for  medium  gallium  concentrations,  up  to  v  1.8 
atomic  percent,  id  depicted  in  figure  2.  From  the  upper  half  of  the 
figure  it  is  apparent  that  considerable  gallium  redistribution  occurs 
during  regrowth.  Redistribution  is  far  more  extensive  than  that  observed 
for  lower  gallium  concentrations  and  the  final  gallium  distribution  has  a 
broad  maximum  distributed  over  v  20  nm  at  the  silicon  surface.  Re¬ 
distribution  is  again  observed  to  occur  at  the  recrystallisation  interface 
but  the  redistribution  peak  broadened  rapidly  towards  the  surface  during 
the  final  anneal.  Impurity  redistribution  of  this  kind  is  generally 
indicative  of  polycrystallite  grain  boundaries  and  hence  this  data  suggests 
that  the  surface  region  of  the  silicon  sample  is  no  longer  amorphous  but 
polycrystalline.  This  is  supported  by  X-ray  texture  camera  analysis  as 
discussed  later.  During  regrowth  the  initial  peak  gallium  concentration 
of  v  1.8  atomic  percent  is  reduced  to  v  1.0  atomic  percent,  a  44 Z  reduction, 
which  is  considerably  more  than  the  23Z  reduction  observed  for  the  lower 
dose  sample.  Following  recrystallisation  a  large  fraction  of  the  implanted 
gallium  is  again  observed  to  reside  on  substitional  lattice  sites.  The 
peak  substitutional  concentration  of  3.0  x  10  Ga/cm  is  7.5  times  the 
maximum  equilibrium  solubility  of  gallium  in  silicon,  a  value  which  is 
comparable  to  that  achieved  by  liquid  phase  epitaxy  during  pulsed  laser 
annealing^,  v  4  x  10^  Ga/cm^.  Comparison  of  figures  1  and  2  shows  that 
redistribution  commences  at  a  similar  gallium  concentration  in  both  cases. 

From  the  lower  portion  of  figure  2  the  average  regrowth  velocity  is  again 
observed  to  be  significantly  enhanced  by  the  presence  of  gallium.  The 
average  regrowth  velocity  of  0,7  nm/min  is  3.5  times  that  expected  for 
uncontaminated  layers  and  the  maximum  velocity  of  v  1  nm/min  is  a  factor 
of  5  times  that  of  uncontaminated  layers.  As  with  the  low  dose  implant 
epitaxial  regrowth  proceeds  at  near  the  undoped  level  for  gallium 


concentration  below  ^  0.2  atomic  percent  and  then  increases  with  increasing 


gallium  concentration  for  concentrations  above  this  value.  It  is  also 
apparent  from  the  regrowth  behaviour  depicted  that  epitaxial  re¬ 
crystallisation  does  not  proceed  to  completion  for  this  sample  but  ceases 
^  20  nm  from  the  surface.  This  is  consistent  with  the  formation  of  a 
polycrystalline  surface  layer  as  suggested  above.  Similar  behaviour 
has  previously  been  reported  for  indium  implanted  (100)  silicon  during 
thermal  annealing  at  525°C^^‘  The  conditions  under  which  a  polycrystalline 
layer  will  form  are  difficult  to  predict  since  the  nucleation  of  such  layers 
is  a  function  of  the  regrowth  velocity,  the  substrate  temperature  and  the 
type  and  concentration  of  impurity,  and  many  of  these  parameters  -are^—  - — 
interelated. 
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The  anneal  behaviour  of  a  high  fluence  gallium  implant  7.5  x  10  Ga/cm  ,  is 
depicted  in  figure  3.  Comparison  of  the  as-implanted  and  post  anneal  gallium 
distributions  shows  that  significant  bidirectional  redistribution  of  gallium 
has  occurred  during  the  total  ancealtime  of  45  minutes  at  465°C.  The  form 
and  the  extent  of  the  gallium  redistribution  is  inconsistent  with  an  epitaxial 
regrowth  process  and  it  again  suggests  the  presence  of  polycrystallite  grain 
boundaries.  The  extent  and  the  bidirectional  nature  of  the  redistribution 
suggest  that  polycrystallites  exist  throughout  the  region  of  the  originally 
amorphous  layer.  It  should  be  noted  that  little  or  no  gallium  redistribution 
was  observed  during  the  first  15  minutes  of  the  anneal  but  after  this  initial 
period  redistribution  occurred  rapidly.  Fletcher  et  al^  have  also  noted 
the  formation  of  an  extensive  polycrystalline  layer  for  high  indium 
concentrations  in  (100)  silicon  during  annealing  at  500°C. 


The. regrowth  behaviour  depicted  in  the  lower  half  of  figure  3  confirms 
that  epitaxy  is  far  from  complete.  The  regrowth  rate  is  observed  to 
be  severely  retarded  ( y  0.15  nm/min)  by  the  presence  of  high  gallium 
concentration  and  comparison  with  figure  2  suggest  that  interfacial 
gallium  concentrations  in  the  range  1-2  atomic  percent  are  sufficient 
to  allow  polycrystallite  nucleation.  Whether  the  polycrystallite 
nucleation  is  a  consequence  of  retarded  epitaxy  or  of  gallium  concentration 
cannot  be  deduced  from  this  data. 

Figure  4  summarises  the  extent  of  epitaxial  recovery  for  the  low,  medium 
and  high  dose  implants  of  figures  1,2  and  3.  The  upper  half  of  the  figure 
shows  Rutherford  backscattering  spectra  for  samples  as-implanted  and  after 
anneal.  The  extent  of  epitaxial  recovery  is  plotted  in  the  lower  half  of 
the  figure  along  with  the  fraction  of  implanted  gallium  residing  on 
substitutional  lattice  sites  after  annealing. 

Following  the  final  anneal  stage  each  annealed  sample  was  examined  in  a 
cylindrical  X-ray  texture  camera  to  establish  the  nature  of  the  annealed 
surface  layer.  X-ray  analysis  was  performed  on  samples  as  removed  from 
the  furnace  and  after  removing  the  thin  native  oxide  5  nm)  which  had 
developed  throughout  analysis.  The  oxide  was  observed  to  contribute  to 
the  diffraction  pattern  and  to  avoid  confusion  only  the  results  for  etched 
samples  are  presented  here.  X-ray  texture  camera  results  are  reproduced 
in  figures  5  and  6,  Figure  5a  shows  the  diffraction  pattern  obtained  from 
crystalline  silicon  for  comparison  and  figures  5b,  6a  and  6b,  show  the 
results  obtained  for  the  low,  medium  and  high  dose  samples  discussed  above. 
For  the  camera  geometry  employed  in  this  study  ’  *  the  vertical  axis  of 

the  photographs  as  displayed,  represents  variable  2  8  ,  where  6,.  is  the 
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Bragg  angle  of  the  reflection.  Horizontal  bands  of  uniform  intensity 
thus  represents  scattering  from  randomly  oriented  planes  or  polycrystals. 


Horizontal  lines  of  variable  intensity  imply  polycrystallites  of  preferred 
crystallographic  orientation. 

Comparison  of  figures  5a,  5b,  6a  and  6b  clearly  suggests  that  preferentially 
oriented  polycrystallites  exist  in  each  of  the  implanted  and  post  annealed 
samples.  The  intensity  of  the  polycrystal  bands,  although  not  accurately 
reproduced  here,  is  observed  to  increase  as  the  extent  of  epitaxial  recovery, 
as  shown  in  figure  4,  decreases.  Thus  the  increasing  intensity  probably 
reflects  the  extent  of  the  polycrystal  layer.  It  is  interesting  to  note 
that  the  orientation  of  the  crystallites  is  the  same  in  all  cases  independent 
of  the  extent  of  the  polycrystal  layer.  This,  along  with  the  fact  that 
a  preferred  orientation  exists,  suggests  that  the  polycrystalline  layer  is 
nucleated  from  the  amorphous-crystal  interface. 

A  particularly  surprising  outcome  of  the  X7ray  analysis  is  the  existence 
of  preferentially  oriented  polycrystals  in  the  lowest  dose  sample.  The 
RBS  spectrum  of  this  sample  shown  in  figure  4,  suggests  complete  epitaxial 
recover,  highlighting  the  insensitivity  of  the  RBS/channelling  technique 
to  certain  crystal  defects, 

SUMMARY  AND  CONCLUSIONS 

The  Rutherford  backscattering/ channelling  technique  has  been  combined  with 
X-ray  texture  camera  analysis  to  characterise  the  recrystallisation  behaviour 
of  gallium  implanted  (100)  silicon  during  furnace  annealing  at  465 °C. 

For  peak  gallium  concentrations  £0.3  atomic  percent  the  RBS/channelling 
technique  suggests  that  epitaxial  recrystallisation  proceeds  to  completion 
with  a  substantial  fraction  of  the  implanted  gallium  being  incorporated  onto 
substitutional  lattice  sites.  The  epitaxial  recrystallisation  rate  is 
enhanced  by  the  presence  of  these  low  gallium  concentrations  and  is  observed 


to  increase  with  increasing  concentration  for  concentrations  £0.2  atomic 
percent.  During  epitaxial  recovery  a  small  fraction  (13%)  of  gallium  is 
redistributed  by  the  advancing  amorphous-crystal  interface.  The  gallium 
incorporated  into  the  recrystallised  layer  is  highly  substitutional  with  a 
peak  substitutional  concentration  2.5  times  the  maximum  equilibrium 
solubility  limit.  RBS/channelling  analysis  suggests  that  epitaxy  is  complete 
but  X-ray  texture  camera  analysis  implies  the  existence  of  a  .very  thin 
preferentially  oriented  polycrystalline  layer.  Comparison  of  the  RBS 
spectra  and  the  X-ray  data  further  suggests  that  this  layer  is  <5  rm  thick. 

When  the  peak  gallium  concentration  is  increased  to  ''•l.  8  atonic  percent 
epitaxial  recrystallisation  no  longer  proceeds  to  the  near  surface  region 
but  ceases  some  20  nm  from  the  silicon  surface.  In  the  region  where 
epitaxial  recrystallisation  occurs  the  regrowth  rate  is  again  observed 
to  be  enhanced  by  the  presence  of  gallium  concentrations  £0.2  atomic  percent, 
despite  the  fact  that  much  of  the  implanted  gallium  is  redistributed  by 
the  amorphous -crystal  interface.  Gallium  incorporated  into  the  regrown 
silicon  layer  is  highly  substitutional  and  a  peak  substitutional  concen¬ 
tration  7.5  times  the  maximum  equilibrium  solubility  limit  was  observed, 
a  value  which  is  comparable  to  that  observed  following  pulsed  laser 
annealing.  Gallium  redistribution  within  the  20  nm  surface  region  suggests 
the  existence  of  polycrystal  grain  boundaries.  X-ray  texture  camera  analysis 
confirms  this  and  further  shows  that  the  polycrystallites  are  preferentially 
oriented. 

15  2 

When  the  fluence  of  40  keV  gallium  is  increased  to  7.5  x  10  Ga/cm  ,  a  peak 
gallium  concentration  of  "M.O  atomic  percent,  annealing  at  465°C  results 
in  little  epitaxial  recrystallisation.  Epitaxy  proceeds  for  ^5  nm  during 
which  the  regrowth  rate  is  severely  retarded.  The  gallium  distribution  is 


observed  to  remain  approximately  Gaussian  during  the  first  15  minutes  at 


465°  and  then  it  rapidly  redistributes  both  towards  the  silicon  surface  and 
towards  the  original  amorphous- crystal  interface.  This  rapid  redistribution 
is  acain  consistent  with  grain  boundary  diffusion  and  X-ray  analysis  again 
confirms  the  existence  of  preferentially  oriented  crystallites. 
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FIGURE  CAPTIONS 


Figure  1: 


Figure  2: 


Figure  3: 


Figure  4: 


Figure  5: 


Figure  6: 


(a)  Gallium  distribution  as  obtained  from  RSS  analysis 
of  as-implanted  (full  curve)  and  post  annealed 
gallium  implanted  (100)  silicon  samples. 


(b)  Amorphous  layer  recovery  as  a  function  of  anneal 
time  at  465°C.  Dotted  line  shows  regrowtji  rate  of 
undoped  amorphous  silicon.  ^ 


Gallium  distribution  and  regrowth  kinetics  as  described 
in  Figure  1. 

7'^  IS  2 

Implant:  40  keV,  -aJ-rJT  x  10  Ga/ca  . 


Gallium  distribution  and  regrowth  kinetics  as  described 


in  Figure  1 
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Imolant: 


b  X  ^  2 

40  keV,  J~c*rx  10^  Ga/cm  . 


(a) 


Channelled  RBS  spectra  as-implanted  and  following 
final  anneal  for  6  x  10“  Ga/cm2  -  0,  3.2  x  lo-1"3  Ga 
cm2  -A,  and  7.5  x  lO1^  Ga/cm2  -□  implants. 


(b)  Amorphous  layer  recovery  and  substitutional  gallium 
fraction  as  a  function  of  gallium  fluence. 


X-ray  texture  camera  diffraction  patterns  for  pure 
crystalline  silicon  (top)  and  for  the  sample  examined  in 
Figure  1  after  final  anneal  (bottom)  . 


X-ray  texture  camera  diffraction  patterns  for  the  sample 
examined  in  Figure  2  (top)  and  Figure  3  (bottom)  after 
final  anneal. 
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ION  IMPLANTATION  DAMAGE  IN  InP 


ABSTRACT 

The  disorder  generated  by  40  keV  light  (N+)  and  heavy  (Bi+)  ion  irradiation 
of  InP  at  40K  and  room  temperature  has  been  measured,  using  Rutherford 
backscattering  channelling  techniques,  as  a  function  of  ion  flux  density  and 
fluence.  For  the  light  ion  irradiation  the  damage  retained  in  the  substrate 
is  highly  dependent  upon  irradiation  temperature  and  upon  flux  density  for 
room  temperature  irradiation.  Such  dependencies  are  much  weaker  for  heavy 
ion  irradiation.  These  results,  .together  with  the  fluence  dependence  of 
disorder,  are  consistent  with  a  mainly  direct  impact  amorphisation,  stable 
against  annealing,  process  with  heavy  ion  implantation  and  a  mainly  simpler 
defect  generation,  unstable  against  annealing,  process  with  light  ion 
implantat i :r . 

Studies  of  disorder  generation  in  both  the  In  and  P  sublattices  are  also 
di scussed . 


^Institute  of  Low  Energy  Nuclear  Physics,  Beijing  Normal  University,  Beijing, 
China. 
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INTRODUaigN 

The  use  of  compound  semiconductors  in  the  electronics  industry  has  grown 
considerably  in  the  last  few  years,  due,  mainly,  to  the  wide  variety  of 
electronic  and  electro-optical  properties  exhibited  by  these  materials. 

Doping  these  compounds  by  conventional  diffusion  techniques  is,  however, 
quite  difficult  since  many  of  the  useful  substrates  have  rather  low  melting 
points  and  the  solubility  of  potentially  promising  dopants  is  often  very 
low.  Ion  implantation  is,  therefore,  a  very  attractive  alternative  for 
doping  compound  semiconductors  and,  indeed,  it  has  already  been  extensively 
used  for  such  applications^. 

Indium  phosphide  and  its  ternary  and  quaternary  compounds  are  currently 

receiving  attention  for  electro-optical  and  electronic  applications,  and 

ion  implantation  is  commonly  employed  as  a  doping  technique  for  these 

materials.  As  with  other  semiconductor  materials,  both  elemental  and 

compound,  ion  implantation  of  InP  can  result  in  severe  radiation  damage.  I 

This  damage  is  normally  detrimental  for  doping  applications,  but  it  is 

extremely  useful  to  know  the  form  of  the  disorder  (i.e.  point  defects, 

extended  defects  or  amorphous  layer),  as  this  will  determine  the  optimum 

(2) 

annealing  procedure  for  damage  removal'  ' .  In  addition,  twere  are  applications 

where  radiation  damage  is  a  desirable  result  of  implantation^,  (e.g.  for 

providing  electrically  isolated  regions  by  compensating  carriers),  and  J 

I 

again  it  is  desirable  to  know  the  form  and  distribution  of  the  disorder. 

Despite  the  relatively  recent  interest  in  InP,  considerable  experimental 
data  exists  for  implantation  into  this  material  Only  recently, 

however,  have  somewhat  exploratory  measurements  of  damage  profiles  and  their 
fluence  dependence  been  reported ^ .  In  such  investigations  only  the  sum 
of  the  In  and  P  displacements  have  been  observed.  In  the  present  study  we 
report  on  the  fluence  and  flux  dependence  of  disorder  in  InP  for  two  cases 


where  quite  different  damaging  mechanisms  may  be  expected  to  operate,  i.e. 
low  and  high  mass  ions.  Investigations  were  carried  out  at  both  room 
temperature  and  40K  and  the  separate  depth  distributions  and  total  displaced 
atom  densities  of  In  and  P  were  observed. 


EXPERIMENTAL 

Commercially  prepared  (100)  InP  samples  were  bombarded  with  40  keV  NT  and  Bi+ 

16  2 

ions  at  room  temperature,  to  doses  up  to  5  x  10  ions/cm  .  A  variety  of 

-2  -2 

current  densities  between  1  ua  cm  and  10  pa  cm  were  employed  for  implant¬ 
ation  so  that  dynamic  annealing  processes  could  be  investigated.  The  target 

_7 

chamber  pressure  was  maintained  at  ^10  torr  during  implantation  and  samples 
were  tilted  ^7°  from  normal  incidence  to  minimise  ion  channelling  effects. 
After  implantation  samples  were  analysed  by  Rutherford  backscattering  and 
channelling  of  2  MeV  He+  ions  on  a  separate  accelerator  at  Salford.  The  ion 
current  densities  indicated  above  are  average  values  for  the  swept  line 
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focus  ion  beam.  Instantaneous  current  densities  are  up  to  a  factor  of  10 
larger  than  these  values. 

In  compatible  experiments  conducted  at  McMaster  University  40  keV  N+  and  Bi+ 
were  implanted  into  InP  at  40K,  at  much  lower  current  densities.  These 
implanted  samples  v/ere  then  analysed,  in  situ,  with  MeV  He+  .ion  backscattering/ 
channelling,  with  the  analysis  beam  derived  from  a  separate  on-line 
accelerator.  This  facility  has  been  described  elsewhere^). 

Backscattering/channel 1 ing  spectra  were  recorded  following  each  implant  and 
random  and  aligned  (100)  axis  data  acquired  for  each  implant  condition. 

It  was  clearly  observed  that  the  analysing  He+  probe  beam  caused  substantial 


disorder  production  in  the  InP,  even  at  room  temperature.  Si  shows  a  similar 
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behaviour  at  low  temperatures^.  Consequently  the  disorder-He+  fluence 
behaviour  was  determined  for  a  non  heavy  ion  implanted  InP  substrate  over 
a  fluence  range  much  larger  than,  but  including,  the  low  fluence  He+  probe 
level  employed  for  heavy  ion  induced  disorder  analysis.  In  this  way  disorder 
introduced  during  analysis  could  be  subtracted  (linearly^)  from  the 
observed  disorder  data.  In  order  to  convert  measured  backscattering  data 
to  displaced  atom  densities  and  to  fit  a  depth  scale  to  the  backscattered 
energy  scale,  standard  routines  of  energy  variable  scattering  cross-section, 
linear  dechannelling  approximation  and  stopping  powers  (deduced  from  tables) 
were  employed. 


?§§yLILAND_DISCUSSION 

Figure  1  displays  typical  RBS/channel  1 ing  data  following  increasing  fluence 

40  keV  N+  implantation  at  room  temperature  and  for  a  N+  flux  density 

_2 

equivalent  to  1  ua  cm  .  The  results  displayed  are  in  fact  smoothed  with 
respect  to  the  raw  data  by  five  point  averaging.  Even  with  this  procedure, 
it  is  clear  that  retrieval  of  displaced  atom  densities  as  a  function  of 
depth  (from  the  increased  backscattering/channelled  yield  as  a  function  of 
channel  number  or  backscatter  energy)  for  P  atoms,  in  particular,  is  difficult 
because  of  the  high  dechannelling  continuum  in  the  P  region  of  the  backscatter¬ 
ing  spectrum.  For  this  reason,  we  have  not  attempted  here  to  derive 
accurate  displaced  atom  densities  for  comparison  between  the  P  and  In  depth 
distributions.  Total,  depth  integrated,  disorder  densities  have  been 
derived  from  data  such  as  in  Figure  1  however,  with  an  estimated  error  in 
this  density,  N^,  for  P  of  +10%.  The  results  of  such  evaluations  are  shown 

in  Figure  2a  where  log  N^,  for  P  and  In  is  plotted  as  a  function  of  log 

-2 

(fluence,  $)  for  a  1  ua  cm  ion  flux.  In  deducing  these  curves  subtractions 


of  the  disorder  produced  by  the  2  MeV  analysing  He+  beam^  ,•  which  was 
observed  to  be  substantial  even  at  room  temperature,  from  the  apparent 
measured  N+  induced  disorder  were  effected.  Within  the  accuracy  of  the 
measurements  Figure  2a  suggests  that  P  and  In  displaced  atom  densities  are 
approximately  equal.  The  less  reliable  depth  distribution  data  also 
suggest  identity  of  the  P  and  In  displaced  atom  densities  at  all  depths. 

Disorder-fluence  data  similar  to  that  of  Figure  2a  was  obtained  at  different 

_2 

ion  flux  densities  of  1  pa,  5  pa  and  10  pa  cm  with  results  displayed  in 
Figure  2b  for  the  more  accurately  determined  In  displacements.  It  is  clear 
that  except  at  the  highest  disorder  levels,  near  quasi-saturation  where 
disorder  increases  only  slowly  with  fluence,  the  effect  of  increasing  ion 
flux  density  is  to  increase  the  disorder.  The  effect  is  equally  pronounced 
in  the  depth  distributions  of  displaced  atom  density  which  are  shown  in 
Figures  3a  and  3b  for  In  displacements  for  N+  fluences  of  5  x  10^  m~2  and 

If  O  O 

10  cm  respectively  for  ion  flux  densities  of  1  pa,  5  pa  and  10  pa  cm 
for  each  fluence.  These  results  clearly  indicate  that  increasing  ion  flux 
density  leads  to  some  increase  in  disorder  at  all  depths  and  that  the 
enhancement  becomes  more  pronounced  at  lower  disorder  levels,  i.e.  the 
effect  is  more  pronounced  at  all  .depths  for  the  lower  fluence  data  of  Figure 
3a  and  more  pronounced  deeper  in  the  solid  in  both  Figures  3a  and  3b. 

Similar  results  to  the  above  were  obtained  using  the  heavier  Bi  ions  at 
40  keV.  Since  disorder  is  concentrated  nearer  the  surface  the  P  displaced 
atom  density  is  even  more  difficult  to  deduce  and  so  results  only  for  In 
will  be  given.  Figure  4  displays  log  N^/log  <$>  data  for  40  keV  implants  at 

room  temperature,  at  low  fluence  levels,  with  an  0.01  pa  cm  flux  density, 
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but  at  higher  fluence  levels  at  1  pa  cm  .  This  strategy  was  required  to 
obtain  accurate  fluence  values  (longer  implant  times)  at  the  lower  fluence 
levels  which  are  demanded  because  of  the  very  rapid  initial  disordering 
rate  (compare  the  high  disorder  level  achieved  with  low  fluence  Bi+  implants 


of  Figure  4  v/ith  the  much  lower  disorder  levels  for  corresponding  fluences 

of  N+  in  Figure  2a)  but  also  to  allow  acceptable  implant  times  at  higher 

fluence  levels.  Figure  4  also  shows  data  obtained  with  40  keV  Bi+  implants 

at  40  K  using  an  on-line  implanter  -  He+  accelerator  analysis  system  at 

McMaster  University  and  described  el sewhere^ .  This  figure  indicates 

that  at  high  fluence,  and  near  quasi-saturati on  disorder,  the  total  disorder 

is  rather  independent  of  substrate  temperature  and  it  was  also  observed 

that  at  high  fluence  (disorder)  levels,  the  disorder  was  almost  flux  density 

independent.  This  behaviour  is  illustrated  in  Figure  5b  which  shov/s  the 

(approximate)  depth  distribution  of  In  displaced  atoms  for  an  implant 
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fluence,  at  room  temperature,  of  5  x  10  cm  and  for  flux  densities  of 

_? 

1  ua,  5  ua  and  10  ua  cm  .  Although  there  is  an  increase  in  disorder  in 
the  deeper  regions  of  the  crystal  with  increasing  flux  density  the  effect 
is  small  but  at  lov/er  fluence  and  disorder  levels  the  effect  is  larger  but 
less  marked  than  for  N+  implants.  Figure  5a  illustrates  this  comparatively 

larger  increase  in  disorder  at  all  depths  with  increasing  flux  density  for 
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a  fluence  of  10  Bi  cm  . 

In  comparing  and  discussing  the  preceding  data  we  note  first  that  the  initial 
rate  of  disorder  production  (from  the  low  fluence,  displaced  atom  intercept) 
is  much  greater  for  Bi+  ions  at  both  room  temperature  and  low  temperature 
than  for  N+  ions.  This  behaviour  is  entirely  consistent  with  earlier 
observations  with  Si  and  GaAs  substrates  and  is  a  reflection  of 

the  high  energy  density  deposition  conditions  induced  by  the  heavier  ion 
which  leads  to  displaced  atom  densities  much  in  excess  of  linear  cascade 
predictions  and  is  the  result  of  largely  direct  impact  amorphisation 
processes^).  The  initial  linear  slope  of  the  log  N^/log  4>  plot  for  low 
temperature  implants  is  a  further  reflection  of  the  probable  dominance  of 
direct  amorphisation  in  individual  cascades  or  spikes^’^'^’^ .  If  the 
present  40  keV  Bi  +  ,  low -temperature  data  is  compared  with  similar  measurements 
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for  GaAs  substrates (similar  mean  mass  target)  it  is  found  that  the 
disorder  at  all  fluence  levels  including  quasi-saturation,  is  approximately 
equal.  It  is  easy  to  understand  the  reason  for  this  similarity  since  cascade 
volumes  and  energy  densities  would  be  expected  to  be  similar  for  the  two 
substrates  whilst  damage  depth  profiles,  which  dictate  amorphous  layer  thick¬ 
nesses^’^  would  be  expected  to  be  similar  also. 
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The  much  lower  disorder  production  rate  for  N+  implantation  suggests,  as  for 
Si  ^»-)  an(j  Qa^s  (10)  substrates,  that  the  major  disordering  process  is  one 
of  simpler  defect  production  with  increasing  defect  accumulation  until  local 
collapse  to  amorphousness  occurs  at  a  critical  defect  density^).  As  these 
amorphous  zones  accumulate  and  overlap  a  continuous  amorphous  layer  results 
as  indicated  by  the  backscattering  data  of  Figure  1.  The  lov/er  displaced  atom 
density  at  the  quasi-saturation  level  for  Bi+  implants  as  compared  to  N+ 
implants  is  a  direct  reflection  of  the  shallower  implant  range  and  disorder 
production  depth  for  the  heavier  ion  which  leads  to  thinner  amorphous  layer 
production. 

It  is  also  important  to  note  that,  to  the  limit  of  present  accuracy,  there 
are  no  major  differences  (for  N+  implantation)  of  either  the  depth  distri¬ 
butions  at  total  displaced  atom  densities  of  the  In  4  P  components.  We  did 
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not  observe,  therefore,  the  theoretically  predicted'  ',  non-stoichiometric 
spatial  distributions  of  In  &  P  vacancies  and  interstitials.  It  should  be 
remembered  that  these  predictions  indicate  only  small  non-stoichiometries, 
within  the  error  limits  of  the  present  data  particularly  for  P,  whilst  the 
predictions  refer  to  different  implant  conditions  (more  energetic  ions  than 
employed  here). 

Turning  next  to  the  effects  of  variation  of  ion  flux  density  it  is  immediately 
apparent  that  for  both  implant  species  and  all  fluence  conditions,  the 
disorder  increases  with-increasing  flux  density.  The  variation  of  disorder 
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with  flux  density  is  quite  generally  an  indication'  '  of  the  operation  of 

defect  migration  and  annealing  processes  which  compete  with  the  generation 

process.  Thermal  annealing  rates  are  expected  to  increase  with  increasing 

substrate  temperature  and  reduce  disorder  levels  and  consequently  the 

increase  in  disorder  with  increasing  flux  density  is  a  clear  indication 

that  the  increasing  flux  density  (and  thus  beam  power  dissipated  in  the 

substrate)  does  not  result  in  sufficient  temperature  increase  to  accelerate 

annealing.  The  process  is  thus  one  v/here  the  increasing  flux  density  more 

successfully  competes  with  thermal  (flux  density  independent)  annealing  to 
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produce  more  stable  disorder'  ' . 

It  is  notable  that  the  effect  of  increasing  flux  density  is  most  marked  for 
the  lighter  ion  implant,  at  the  lower  fluence  and  disorder  levels  for  both 
ion  species  and  in  the  deeper  components  of  the  disorder-depth  distributions. 

i 

In  all  these  conditions  the  nature  of  the  disorder  is  expected  to  be  of 
simpler  defect  form,  even  in  the  case  of  the  heavy  ion  irradiation  where, 
in  the  deep  tails  of  individual  collision  cascades  the  energy  deposition 
density  will  be  insufficient  to  generate  spike  conditions  and  at  the  radial 
peripheries  of  such  individual  cascades  nearer  the  surface  the  damage  will 
be  simpler  before  complete  amorphous  layers  have  accumulated. 

It  is  thus  reasonable  to  conclude  that  the  flux  density  influence  on  measured 
disorder  production  is  a  clear  indication  of  the  production  and  subsequent 
migration,  at  room  temperature,  of  simple  defects.  In  earlier  studies  of 
disorder  production  in  InP,  under  mainly  low  energy  deposition  density 
conditions,  there  are  also  signs  of  annealing  processes,  both  from  the  post¬ 
irradiation  long  term  annealing  at  room  temperature  studies^  and  the  rapid 

(2  3  5 ) 

increase  with  temperature'  *  ’  '  near  room  temperature  of  the  fluence  required 
for  layer  amorphisation.  Annealing  clearly  occurs  under  all  irradiation 
conditions  employed  here,  to  a  greater  or  less  extent,  but  is  most  important 


of  simpler  defect  production  mechanisms  for  lower  energy  deposition  density 
(and  associated  local  defect  density)  conditions  and  direct  impact  amorphisa- 
tion  (or  some  similar  form  of  thermally  stable  disorder)  for  higher  energy 
deposition  density  conditions.  The  existence  of  some  annealable  disorder, 
even  with  the  heavy  ion  (higher  energy  density  deposition)  is  not  only  evident 
in  the  dependence  of  disorder  on  flux  density  shown  in  Figure  5  but  also  in 
the  increase  in  the  low  fluence  disorder  with  decreasing  substrate  temperature 
shown  in  Figure  4.  Similar  low  temperature  studies  for  N+  implants  are 
currently  in  progress  and  will  be  reported  elsewhere  and,  in  view  of  the  very 
significant  influence  of  flux  density  for  this  ion  species,  are  expected  to 
reveal  very  large  increases  in  disorder  generation  rate  when  compared  to  the 
room  temperature  data  of  Figure  2. 

In  view  of  the  observed  flux  density  effects  and  the  implied  existence  of 
annealing  processes  it  is  clear  that  measurements  of  disorder  generation 
rates  per  ion  (N^*)  should  only  be  attempted  with  low  temperature  substrates,  a 
reported  for  Si  ^>9)  an(j  q3^s  (10)  -jn  orcjer  to  compare  with  model  (e.g. 
linear  cascade  or  spike)  predictions.  It  is  for  this  reason  that  detailed 
quantitative  evaluations  of  N^*  have  not  been  attempted  here  but  only 
qualitative  comparisons  between  Bi+  and  N+  implants. 


CONCLUSION 

This  study  has  revealed  that  probable  disordering  mechanisms  resulting  from 
both  light  and  heavy,  relatively  low  energy,  ion  impact  on  InP  are  similar 
to  and  fit  into  the  general  scheme  proposed  for  both  Si  and  GaAs^^. 

Thus  implantation  resulting  in  high  energy  deposition  density  conditions 
leads  predominantly  to  direct  impact  amorphisation  processes  whereas  lower 
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accumulates  and  transforms  to  amorphousness  at  critical  levels.  These  hypo¬ 
theses  are  supported  by  the  strong  influence  of  ion  flux  density  upon  stable 
disorder  production  for  lighter  ion  implantation  and  the  weaker  but  evidently 
similar  effects  with  heavy  ion  implantation,  and  in  this  latter  case  the 
influence  of  substrate  temperature. 

It  is  clear  that  in  technological  applications  careful  control  of  ion  flux 
density  and  substrate  temperature  should  be  exercised  and  that  if  amorphous 
layer  structures  are  desirable  high  flux,  low  temperature  conditions  are 
advised  but  if  such  structures  are  undesirable  low  flux,  elevated  temperature 
conditions  are  advised. 

Within  the  limits  of  experimental  accuracy,  no  major  differences  in  the 
stoichiometry  of  In  and  P  displacements  were  observed. 
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FIGURE  CAPTIONS 


Fig.  1  Random  and  aligned  RBS  spectra  of  2  MeV  He+  ions  from  40  keV 

N+  implanted  InP  at  room  temperature  for  different  incident 
N+  fluences. 

t 

Fig.  2a  The  number  of  displaced  In  &  F  atoms  (N^)  in  InP  as  a  function 

of  40  keV  N+  ion  implantation  fluence  at  room  temperature. 

_2 

The  ion  current  density  was  equivalent  to  1  pa  cm  . 

Fig.  2b  The  number  of  displaced  In  atoms  in  InP  as  a  function  of  40  keV 

N+  ion  implantation  fluence  at  room  temperature.  Ion  current 

-2 

densities  equivalent  to  1 ,  5  and  10  ua  cm  were  employed. 

Fig.  3a  The  effective  depth  (channel  no.)  distribution  of  displaced 

In  atoms  in  InP  following  40  keV  N+  ion  implantation  at  room 
temperature  to  a  constant  fluence  of  5  x  1014  ions  cm-2  at 
equivalent  flux  densities  of  1 ,  5  and  10  pa  cm  . 

Fig.  3b  The  effective  depth  (channel  no.)  distribution  of  displaced 

In  atoms  in  InP  following  40  keV  N+  ion  implantation  at  room 

16  -2 

temperature  to  a  constant  fluence  of  10  ions  cm’  at 

_2 

equivalent  flux  densities  of  1 ,  5  and  10  pa  cm  . 

Fig.  4  The  number  of  displaced  In  atoms  in  InP  as  a  function  of  40  keV 

Bi+  ion  implantation  fluence  at  room  temperature  and  at  40  K. 

Fig.  5a  The  effective  depth  (channel  no.)  distribution  of  displaced  In 

atoms  in  InP  following  40  keV  Bi+  ion  implantation  at  room 

temperature  to  a  constant  fluence  of  10^4  ions  cm’2  at  equivalent 

_2 

flux  densities  of  1,  5  and  10  pa  cm  . 


The  effective  depth  (channel  no.)  distribution  of  displaced 
In  atoms  in  InP  following  40  keV  Bi+  ion  implantation  at  room 
temperature  to  a  constaht  fluence  of  5  x  10^4  ions  cm'2  at 
equivalent  flux  densities  of  1 ,  5  and  10  ya  cm'2. 
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DISORDER  PRODUCTION  AND  ANNEALING  DURING  He  ION  RBS/CHANNELLING  ANALYSIS 
OF  InP 


ABSTRACT 

•  _  _  + 

Studies  of  the  interaction  of  2  MeV  He  ions,  employed  for  Rutherford 
backscattering/channelling  of  disorder  in  InP,  with  the  disorder 
generated  by  pre-implantation  of  40  keV  N+  ions  at  room  temperature 
are  reported.  It  is  shown  that  for  initially  undamaged  substrates, 

He+  irradiation  generates  disorder  whilst  for  N+  implantation  damaged 
substrates  He+  irradiation  anneals  disorder.  Possible  mechanisms  to 
account  for  this  behaviour  are  discussed  and  the  potential  problems 
associated  with  accurate  disorder  measurements  using  RBS/channelling 
outlined. 

INTRODUCTION 

Rutherford  backscattering/channelling  analysis  using  MeV  He  ion  beam  probes 
is  a  common  technique  for  disorder  assessment  in  semiconductors  .  It  is 
known,  however,  that  the  technique  is  perturbational  and  is  generally  unable 
to  assess,  for  this  and  other  fundamental  reasons,  disorder  densities  less 
than  about  1%  total  lattice  randomisation.  Two  of  the  important  perturbational 
processes  are: 

1)  the  generation  of  disorder  by  the  He  probe  additional  to  that  which 
may  already  exist  in  the  substrate,  and 

2)  the  annealing  of  disorder,-  which  may  already  exist  in  the  substrate,  by 
the  He  probe . 

2  3  4  5 

The  first  process  has  been  studied  in  some  detail  recently  for  Si  '  ' 

since  the  effect  is  f&r  stronger  near  the  surface  than  would  be  expected  from 
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the  usual  modified  Kinchin-Pease  damage  generation  model.  Investigations 

with  GaAs  have  also  revealed  the  higher  than  expected  disorder  production 

rate  and,  in  addition,  have  suggested  that  the  disorder  generated  by  the  He 

probe  is  non-linearly  additive  to  that  created  by  a  previous  heavier  ion 

irradiation.  The  second  process  was  noted  in  Si  '  and  GaAs  “  and  more 

recently  annealing^  of  Si+  implanted  Si  on  sapphire  has  been  reported  at 

temperatures  from  200 °C  to  400°C  using  He+,  N+  and  Ne+  ions  in  the  hundreds 

+  14  • 

of  keV  range.  Early  studies  with  Te  implanted  GaP  revealed  no  measurable 
annealing  with  1  MeV  He+  ions  (indeed  the  increased  disordering  process  was 
noted)  nor  with  10  keV  or  15  keV  He+  ions  which  should  possess  defect 


generation  depth  profiles  overlapping  the  Te  disorder  profile. 


As  part  of  a  programme  of  investigation  of  heavy  ion  disorder  production  in 
InP1^  using  He+  ion  backscattering  analysis  we  have  observed  the  operation 
of  both  disordering  and  annealing  processes  with  some  rather  unusual  and 


unexpected  results  which  will  be  presented  and  discussed  in  this  communication 


EXPERIMENTAL 


Full  details  of  the  experimental  technique  have  been  given  elsewhere  and 
only  a  brief  summary,  relevant  to  the  present  study,  will  be  given  here. 
Commercially  prepared  (IOC)  InP  samples  were  bombarded  at  room  temperature 

with  40  keV  N+  ions  to  fluences  of  lO1^  cm  ^ ,  3  x  10^  cm  ^ ,  10^4  cm  ^ , 

15  -2  15  -2  -2 

10  cm  and  5  x  10  cm  at  fixed  dose  rate  of  1  pa  cm  .  This  latter 


parameter  was  controlled  since  it  is  known  to  influence  the  amount  of 
disorder  produced  and  retained  in  the  InP^.  The  samples  were  tilted  'W° 
from  normal  incidence  to  minimise  N+  ion  channelling  effects. 
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After  implantation  the  damage  induced  in  the  InP  was  analysed  by  Rutherford 
scattering/ channelling  studies  using  2  MeV  He+  ions  in  randcxn  or  (lOO)  alignei 
incidence,  168°  scattering  angle  geometry. 

Backscattering  spectra  for  aligned  incidence  were  determined  as  a  function 

of  incident  He+  integrated  ion  charge  up  to  ^50  yC  (equivalent  to  a  He+  ion 

16  —2 

fluence  density  of  <\,3  x  10  cm  )  .  The  minimum  charge  required  for  good 

statistical  evaluation  of  disorder  was  approximately  1.5  yC.  In  subsequent 

evaluation  of  the  measured  data,  the  backscattering  yield  spectra  for  given 

integrated  charge  were  normalised  to  the  equivalent  values  for  a  1.5  yC 

integrated  charge.  He+  ion  currents  of  4.5  nA  were  generally  employed  (over 
2 

an  area  of  M  mm  )  but  some  studies  were  conducted  with  25  nA  ion  currents. 

In  addition  to  use  of  the  He+  ion  probe  for  analysing  N+  implant  induced 
disorder,  studies  were  also  made  of  the  effects  of  He+  probe  beams  on  non- 
implanted  samples  in  which  the  same  strategy  as  above  was  followed  but  with 
an  unimplanted  InP  crystal. 


RESULTS 

Figure  1  shows  backscattering/channelling  spectra  from  non-implanted  InP 
normalised  to  expected  yields  for  1.5  yC  integrated  charge  as  a  function  of 
increasing  He+  ion  charge  (fluence) .  Only  the  portion  of  the  spectrum  in 
the  region  of  backscattering  from  displaced  In  atoms  is  displayed  since 
backscattering  yields  from  the  lighter  mass  P  are  much  lower  and  difficult 
to  determine  accurately  against  the  rising  dechannelling  background  at 
lower  scattering  energies^.  It  is  clear  that,  for  this  initially  undamaged 
substrate,  the  He  probe  ion  generates  In  displacements  at  all  depths,  but 


i n  particular,  near  the  surface  where  a  low  Initial  level  of  In  displacements 
already  exists.  The  increase  in  this  In  disorder  is  determined  from 

the  integrated  area  of  the  peaks  in  Figure  1  with  normal  linear  dechannelling 
corrections  applied16.  An  effective  cross-section  for  disordering  may  be 
deduced  from  .the  initial  slope  of  a  linear  or  semi -logarithmic  plot  of  In 
disorder  as  a  function  of  fluence  from  the  relation  : 


(1) 


dNd 

where  N  is  the  initial  area  of  the  surface  peak  and  — —  is  the  slope  of  the 
o  a® 

increase  in  surface  peak  area  with  respect  to  fluence  plot.  Determination  of 
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this  value  from  the  data  of  Figure!  yields  %  2.5  x  10  cm  ion  .  This 

value  is  a  factor  of  about  four  larger  than  cross-sections  measured  at  40  K 

+  5 

and  300  K  for  1  MeV  He  irradiation  of  Si  ,  but  rather  similar  to  cross-sections 

8  4- 

measured  at  50  K  for  2  MeV  He  irradiation  of  GaAs. 


This  value  of  for  InP  must  however  be  regarded  as  an  upper  limit  since  after 

the  first  4-5vjC  integrated  He+  probe  charge  there  is  little  further  change  in 

the  area  of  the  surface  peak  once  the  dechannelling  correction  has  been  applied. 

Moreover  repetition  of  the  experiment,  but  with  lower  initial,  He+  fluences  showec 

an  initial  negative  value  of  (i.e.  an  annealing  occurred),  with  a  cross- 

section  of  order  10  15  an2.  For  He+  fluences  greater  than  about  5  x  1014  cm  2 

the  process  changed  to  a  disordering  mechanism  with  a  cross-section  similar 

to  the  value  indicated -above.,  but  again  at  still  higher  fluences  (xlO1^  cm  2)  , 

o.  tended  towards  zero, 
d 


The  concurrent  studies  of  heavier  (N  and  Bi+)  ion  implantation  into  InP 
have  shown  that  with  increasing  40  keV  N+  ion  fluences  the  total  In  disorder 
increases  rapidly  in  the  fluence  range  1012  cm  2  to  about  3  x  1014  cm-2  and 
then  much  more  slowly  with  increasing  fluence.  The  RBS/ channelling  data 
reveal  that  in  the  lower  fluence  regime  the  In  peak  does  not  reach  the  random 
level  (i.e.  the  near  surface  region  is  disordered  but  apparently  not  fully 
randomised)  but  at  higher  .fluences  the  In  peak  reaches  the  random  level  and 
then  continues  to  broaden.  In  common  with  earlier  interpretations  for  other 
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amorphised  over  a  small  depth  near  the  surface  and,  as  fluence  is  further 
increased  the  amorphised  layer  thickness  increases.  As  a  result  of  this  data 


‘ 

N+  ion  fluences  were  implanted  to  correspond  to  disorder  levels  from  %lo% 
initial  randomisation,  through  initial  randomisation  and  beyond  this  level  to 
where  the  randomised  layer  thickness  was  increasing  slowly  with  fluence. 

This  required  ion  fluence  densities  to  be  chosen  from  10^  cm  ^  to  5  x  10^ 

-2  +  • 

cm  as  indicated  earlier.  Following  each  of  these  N  implants  the  InP  was 

analysed  with  2  MeV  He+  ions  under  aligned  conditions  and  the  backscattering 
data  determined  for  increasing  He  fluence.  Examples  of  the  results  of  these 
analyses  are  shown  in  Figures  2,  3  and  4  for  initial  N+  ion  implant  fluences 
of  10^,  lO1^  and  5  x  10^  cm  ^  which  correspond  to  situations  of  less  than 
initial  randomisation,  the  region  of  layer  randomisation  and  of  randomised 
layer  thickening  respectively.  These  figures  all  indicate  that,  with  increas¬ 
ing  He  irradiation  the  number  of  displaced  In  atoms  in  the  near  surface  peak 
is  decreased.  For  the  two  lower  N+  implant  fluences  there  is  also  a 
decrease  in  the  behind  the  In  peak  and  a  decrease  in  the  deeper  dechannel¬ 

ling  level.  In  the  case  of  the  highest  fluence  N+  implant  there  is  an  increase 
in  Xmin  and  deeper  dechannelling  accompanying  (as  in  the  lower  N+  fluence 
cases)  a  reduction  in  the  In  peak  height.  In  Figures  2  and  3,  there  is 

clearly  a  narrowing  of  the  In  peak  also  from  the  deeper  side  of  the  peak 

15  -2 

(the  apparent  shift  of  the  surface  side  of  the  peak  for  the  lo  cm  may 
well  be  due  to  a  small  shift  in  the  energy  scale  for  one  particular  analysis 
after  15  yC) ,  whereas  there  is  no  clear  inward  shift  or  peak  narrowing  for 
the  5  x  10  ^  N+  implant  condition.  When,  however,  it  is  remembered  that  the 
dechannelling  increases  with  He+  ion  fluence  in  this  case  then  when  the 
dechannelling  subtraction  is  made  from  the  peak  region,  the  deep  edge  of  the 
peak  does  indeed  move  towards  the  surface  and  the  peak  narrows. 

All  of  the  above  experiments  were  repeated  several  times  with  essentially 
qualitatively  similar  results  although  the  magnitudes  of  the  effects  varied 
slightly.  Very  similar  results  were  obtained  with  He+  ion  probe  currents  of 
4.5  nA  and  25  nA  and  the  results  of  Figure  4  were  also  repeated  when,  following 
data  acquisition  for  15_  uC  He+  collected  charge,  the  counters  were  cleared 


and  an  analysis  performed  for  a  further  1.5  UC  collected  charge.  Dead  time 


6 


magnitude  of  the  effect. 


Figures  2,  3  and  4  all  demonstrate  that  the  effect  of  the  He  ion  probe  is, 
for  all  levels  of  N+  implantation  induced  disorder,  to  induce  an  annealing 
behaviour.  In  the  same  way  as  a  He+  ion  disordering  cross-section  was 
determined,  it  is  possible  to  deduce  an  annealing  cross-section 


where  N,  is  the  initial  level  of  N  induced  disorder, 
a . 

i 


Evaluation  of  a  from  the  data  of  Figures  2,  3  and  4  yields  a  value  of  1-* 
a  *" 

„  -17  2.  -1 

2  x  10  cm  ion  .  This  value  of  a  was  relatively  independent  of  N 

^  d . 

l 

from  ^10%  randomisation  to  above  initial  randomisation  and  did  not  appear  to 

vary  substantially  with  increasing  He+  ion  fluence.  In  the  case  of  N+ 

8  + 

implanted  GaAs  at  40  K  the  effect  of  the  He  probe  was  to  increase  the 

disorder  with  a  cross-section  which  was  relatively  independent  of  . 

i 

In  addition  to  these  measures  of  cross-section,  estimates  may  also  be  made 
of  cross-sections  from  the  increase  or  decrease  of  the  dechannelling  minimum 
immediately  behind  the  surface  peaks  from  a  formalism  equivalent  to 
equations  (1)  and  (2),  i.e. 


(XJ 

d  m 


Deduction  of  these  a  values  from  Figures  1  to  4  yields: 


n(X  )  ^  6  x  10  ^  cm^  ion  ^  for  uniraplanted  InP 
c*  in 


c_(X  )  ^  1.4  x  10  ^  cm^  ion  ^  for  N+  implanted  InP  for  fluences  up  to 
a  m  s c  _9 

about  lO13  cm 
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3  x  10  cm  ion  for  N  implanted  InP  for  a  fluence  of 

5  x  1015  cm-2. 


DISCUSSION 


Considering  first  the  disordering  induced  by  He  irradiation  of  unimplanted 

InP  it  is  noted  that  the  cross-sections  for  Si,  GaAs  and  InP  are  all  of 

3-5 

similar  order.  This,  as  has  been  argued  elsewhere  ,  is  always  substan¬ 
tially  larger  than  would  be  expected  from  calculations  of  damage  production 

7 

cross-sections  based  upon  Kinchin-Pease  calculations  which  account  for 
recoil  multiplication  processes  in  addition  to  direct  He-In  collisions. 
This  phenomenon  is  not  well  understood  but  may  result  from  lower  than  bulk 
value  displacement  energies  for  near  surface  atoms'’  and  from  relaxed  atoms 


around  displacements  contributing  to  He  backscattering 


There  is  some 


evidence  for  the  latter  or  similar  processes  in  the  present  data. 


If  scattering  occurs  only  from  fully  displaced  atoms  then,  as  shown  by  Bf6gh  , 

should  be  linearly  related  to  N^.  This  relationship  has  been  well 
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demonstrated  for  heavier  than  He  ion  implantation  of  Si  ,  GaAs  and  InP 

irradiated  at  low  temperature.  The  present  data  which  indicate  o  ^  2.5 

2.5 

suggest  a  relationship  of  the  form  Xm  *  kN^  ’  .  The  earlier  studies  of  1  MeV 
He+  irradiated  Si^ at  40  K  also  suggest  a  slightly  superlinear  relationship 
between  Xffi  and  whilst  measurements  with  heavier  ion  implanted  Si^  at  40  K 
indicate  a  slight  superlinearity  for  P+  implantation. 


All  of  these  results  suggest  that  strain  effects  are  also  contributing  to 


the  measured  dechannelling  process.  These  may  result  from  atomic  relaxations 


around  completely  displaced  atoms  (and  vacancies)  and  from  defect  clusters 

which  may  form  in  the  present  studies  at  room  temperature  where,  our  parallel 
2^  suggest 

studies  with  heavier  implant  ions'',  considerable  defect  migration  must  occur. 

Turning  now  to  the  He+  probe  annealing  effects  in  already  N+  implant  damaged 
InP  it  is  notable  that,  for  all  disorder  levels,  a  is  considerably  smaller 
than  but  for  the  highest  disorder  level  although  annealing  of  the  peak 
disorder  occurs,  there  is  a  simultaneous  disorder  production  as  evidenced  by 
the  increase  in  with  He+  fluence  (i.e.  is  positive) .  Moreover  whilst 
it  is  true  that  seme  of  the  peak  annealing  occurs  by  apparent  regrowth  of 
the  "amorphous"  layer  from  the  crystalline  substrate  as  evidenced  by  the 
shift  towards  the  surface  of  the  inner  boundary  of  the  peak,  it  is  also  clear 
that  annealing  occurs  within  this  "amorphous"  layer  itself  since  the  peak 
height  is  everywhere  continuously  reduced  by  increasing  He+  irradiation. 

The  mechanisms  responsible  for  this  unusual  behaviour  are  currently  unclear 
but  seem  not  to  be  associated  simply  with  thermal  reordering  within  the 
randomised  layer  and  at  the  random-crystal  interface  since  changes  in  He+ 
beam  power  by  a  factor  of  5  resulted  in  no  notable  differences  in  anneal 
behaviour.  We  speculate,  however,  that  although  the  RBS/channelling  data 
reaches  the  random  level  that  this  is  not  a  precise  indicator  of  full 
amorphisation  and  that  the  near  surface  region  may  be  either  totally  micro¬ 
crystalline  or  largely  amorphous  with  small  included  crystalline  regions. 

If  either  were  the  case  then  some  recrystallisation  in  this  layer  would  be 
possible  by  He+  induced  defect  generation  and  migration  and/or  by  inelastic 
energy  loss  processes  which  disrupt  the  rearranged  bonding  configuration 
created  by  implantation.  We  propose  to  further  investigate  these  processes 

by  variation  of  probe  energy,  and  species  in  order  to  vary  elastic  and 

8 

inelastic  energy  loss  rates  as  has  been  initiated  with  C-aAs.  If  such  local 
reordering  does  in  fact  occur  than  it  may  partly  account  for  the  positive 
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0j(X  )  observed  for  initially  heavily  disordered  layers  since  such  crystallite 
zones  and  their  environs  could  lead  to  local  strains  and  enhancement  of 
dechannelling.  It  is  also  clear  that  microscopic  (TEM)  and  low  incidence 
angle  diffraction  techniques  would  be  helpful  in  obtaining  a  clearer  under-  • 
standing  of  the  nature  of  the  near  surface  damaged  layer. 

9-13  + 

It  is  thus  believed,  in  common  with  earlier  suggestions  ,  that  the  He 
probe  annealing  process  is  due  to  simple  defect  production  which  interacts 
with  the  more  stable  initial  disorder  and  reduces  and  modifies  this  disorder. 

Finally  the  present  results  indicate  a  problem  which  has  not  arisen  previously 

in  the  He+  probe  analysis  of  Si  and  GaAs .  In  these  materials  it  has  been 

noted  that  the  He+  probe  always  creates  additional  disorder  to  that  induced 

8  + 

by  implantation  and  even  if  these  disorders  are  non-linearly  additive  the  He 

probe  effect  can  be  described  by  a  more  or  less  constant  production  cross- 

section  process.  Appropriate  deconvolution  of  measured  disorder  data  to  give 

5  8  * 

actual  disorder  data  can  then  be  made  '  .  In  the  case  of  InP  it  appears  that 

for  disorder  levels  somewhere,  as  yet  poorly  specified,  between  1%  and  10% 

probably 

initial  randomisation,  O’ Changes  sign.  At  high  disorder  levels,  and  since  He 
probe  fluences  necessary  for  analysis  are  generally  low,  the  sign  and  low 
magnitude  of  o  result  in  little  adjustment  to  be  made  to  measured  data  to 
recover  actual  disorder  data.  Such  adjustments  are  much  more  severe  in  the 
low  implant  fluence  (low  disorder  concentration)  regimes  where,  as  just  noted 
a  values  are  poorly  described  in  sign  and  magnitude.  Thus  the  accurate 
evaluation  of  low  disorder  concentrations  in  InP  is  much  poorer  than  either 
Si  and  GaAs.  At  a  fundamental  level  the  variation  of  a  with  increasing 
disorder  can  be  explained  in  terms  of  a  disorder  production  cross-section 
which  decreases  with  and  a  competing  disorder  annealing  cross-section 
o_  which  either  increases  or  remains  constant  with  increasing  N  .  Evidence 

a  Q 

for  some  competitive  effect  exists  from  the  data  of  Figure  4  discussed  earlier 


which  reveals  both  a  oa  component  in  the  disorder  peak  annealing  and  a  ( y^) 

component  in  the  dechannelling  increase.  Again  we  can  only  speculate  on  the 

reason  for  this  behaviour  at  present.  Thus  for  an  unimplanted  substrate  N^ 

may  represent  In  atoms  in  either  or  both  a  relaxed  state  due  to  surface  ’ 

proximity  and  in  a  surface  oxide.  After  implantation  N,  represents  In  atoms 

fli 

displaced  and  partly  relaxed  from  lattice  sites  as  a  result  of  the  implantation 
induced  collision  cascades.  The  interaction  of  the  He  beam  with  these  diff¬ 
erent  configurations  of  In  atoms,  and  indeed  P  atoms,  may  well  be  responsible 
for  the  observations  together  with  processes  which  may  be  non-linear  in  N^ 
such  as  defect  migration  and  annihilation  with  existing  defects  which  will 
increase  in  importance  with  increasing  N^.  It  is  quite  clear  that  further 
studies  of  these  disordering  and  annealing  processes  induced  by  He+  (and  other 
probe)  species  following  a  variety  of  ion  implant  conditions  are  necessary 
in  order  to  optimally  employ  RBS/channelling  as  an  implant  disorder  analysis 
tool  with  InP. 
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channel  number 

Fig.  f  Rutherford  backscattering /channelling  spectra  of  random 
and  <fOO>  aligned  incidence  2MeV  He+  Ions  from  In  atoms 
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Fig.  2  Rutherford  backscattering /channelling  spectra  of  random  and 
<400)  aligned  incidence  2MeV  He+  ions  from  In  atoms  in  an  InP 
substrate  implanted  with  40keV40u  N +  ion&  cm"2  as  a  funr*>w 
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Fig.  3  Rutherford  backscattering/channelling  spectra  of  random  and 
<fOO>aligned  incidence  2MeV  He+ions  from  In  atoms  in  an  InP 
..  _ _ substrate  implanted  with  4QkeV  fQ<5  N*  ions  cm”2  as  a  functio 


General  Conclusions 


Over  the  period  of  the  grant  award  the  investigators  at  Salford 
have  undertaken  two  functions.  Firstly  an  analytical  assessment 
of  samples  provided  by  the  sponsor  at  Wright  Patterson  AFB  have 
been  undertaken  to  the  sponsors  requi rements .  The  results  of 
this  work  were  described  in  our  earlier  report  and  not  only 
reveal  problems  associated  with  strain  at  silicon  on  sapphire 
interface  structures  but  also  suggested  further  studies  which 
could  be  performed  at  Wright  Patterson. 

Secondly  an  associated  research  programme  on  ion  implantation 
processes  in  Si,  GaAs  and  more  recently  InP  has  been  undertaken 
to  provide  basic  intelligence  to  the  potential  application  of 
the  ion  implantation  technique  in  device  construction  relevant 
to  the  interests  of  the  sponsor.  The  results  of  these  studies 
have  been  documented  in  the  earlier  and  present,  final  reports. 
It  is  believed  that  considerable  new  understanding  of  the 
fundamental  processes  of  damage  creation  attendant  upon  ion 
implantation  of  these  semiconductor  materials,  and  upon  impurity 
incorporation  during  furnace  annealing  of  Si,  has  been  achieved. 
This  is  beneficial,  not  only  to  the  grant  sponsor  and  the 
University  investigators,  but  to  the  international  scientific 
and  technological  community  concerned  with  the  use  of  ion 
implantation  in  the  present  and  potential  applications  to 
semiconductor  device  production.  The  University  investigators 
express  gratitude  to  the  USAF  for  provision  of  the  financial 
support  which  has  enabled  prosecution  of  these  investigations. 


